Self-energy feedback and frequency-dependent interactions in the functional 
renormalization group flow for the two-dimensional Hubbard model 
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We study the impact of including the self-energy feedback and frequency-dependent interactions 
on functional renormalization group grows for the two-dimensional Hubbard model on the square 
lattice at weak to moderate coupling strength. Previous studies using the functional renormalization 
group had ignored these two ingredients to large extent, and the question is how much the flows 
to strong coupling analyzed by this method depend on these approximations. Here we include the 
imaginary part of the self-energy on the imaginary axis and the frequency-dependence of the running 
interactions on a frequency mesh of 10 frequencies on the Matsubara axis. We find that i) the critical 
scales for the flows to strong coupling are shifted downwards by a factor that is usually of order one 
but can get larger in specific parameter regions, and ii) that the leading channel in this flow does 
not depend strongly on whether self-energies and frequency-dependence is included or not. We also 
discuss the main features of the self-energies developing during the flows. 

PACS numbers: 71.10.Fd, 71.10.Hf, 74.20.Mn 
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I. INTRODUCTION 



Weakly to moderately correlated, itinerant electron 
systems like iron arsenide superconductors or carbon ma- 
terials are in the center of todays condensed matter re- 
search. Nevertheless, a detailed quantitative description 
of their low-temperature phases is often not available. 
Two potential obstacles are that these materials, de- 
spite not being strongly correlated, can still host var- 
ious competing trends that invalidate simpler theoreti- 
cal approaches, and that the underlying lattice causes 
anisotropic behavior that can lead to significant quanti- 
tative effects such as anisotropic quasiparticle properties 
or anisotropic energy gaps in the ordered states. The 
functional renormalization group (fRG) method has be- 
come a widely used tool for the investigation of weakly 
coupled fermions on low-dimensional lattices (for a re- 
cent review, see^) as it helps in both these aspects. In 
contrast to often employed random-phase approximation 
(RPA) or mean field approaches its is capable of de- 
scribing competing interactions, and it is also capable 
of resolving stronger momentum-dependencies of the ba- 
sic obervables like pairing gaps etc. While fRG applica- 
tions in two dimensions mainly focussed on the Hubbard 
model on the square lattice^—, more recent studies ana- 
lyze also multiband models, in particular for iron arsenide 
superconductors^— and graphene systemsSril. The bulk 
of these fRG works uses an approximation to the full 
flow equations in one which neglects the self-energy com- 
pletely. The main object is interest in the effective, flow- 
ing interaction vertex, whose strongest components give 
information on the leading correlations in the effective 
theory at lower scales. Only in a few works, either the 
flow of the self-energy was computed without allowing 
for a feedback on the flow of the interaction verte x 12 ' 13 , 
or the self-energy was parameterized using a quasiparti- 
cle weight^—. These studies gave the following infor- 



mation: a) The fRG flows with neglected self-energy are 
good approximations in the sense that in typical cases the 
leading low-energy instabilities are not altered in their 
existence and character by including (parts of) the self- 
energy, b) A detailed study of the self-energy suggests the 
occurrence of interesting and observable effects such as 
anisotropic quasiparticle degradation, dispersion renor- 
malizations and partial gap openings. Another approx- 
imation that has been used in most applications of the 
fRG is to neglect the frequency-dependence of the in- 
teraction vertex. This dependence is usually absent in 
the initial condition of the flow, i.e. the bare electron- 
electron interaction is not retarded. But during the flow, 
the subsequent inclusion induces a frequency-dependence 
that can cause at least quantitative effects. Again, the 
frequency-dependence of the vertex has been taken along 
in a handful of studies^—. The upshot here is again that 
the previous flows without frequency-dependence are not 
so bad and not drastic changes of the leading instabilities 
are expected when the approximation is improved. 

All these rather qualitative results are encouraging, 
but in order to develop the fRG into a quantitative 
method that is able to describe materials more directly, 
the approach has to be improved. One group of work- 
ers have embarked on this mission using ansatzes for 
the frequency-dependences of vertex functions and self 
energie s 20 : 21 . Here we want to use a more direct ap- 
proach. We treat the frequency-dependence of the self- 
energy and the interaction vertex by discretizing the Mat- 
subara axis in N u — 10 patches, in a similar spirit as the 
wavevector-dependence is patched into Nk = 32 patches 
around the Fermi surface. This way, the interaction 
vertex of the one-band model becomes an object with 
(Nk x N^) 3 components that can still be dealt with nu- 
merically. Also, the imaginary part of the self-energy is 
computed with some resolution, given by N u values on 
the Matsubara axis. This allows us, e.g., to test how 
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adequate it is to use a simple Z-factor in order to de- 
scribe the flowing self-energy. The self-energy is then fed 
back into the flow of the interaction vertex by using full 
Green's functions on the internal lines of the one-loop 
diagrams in their flow. 

The main questions we can ask using this refined 
scheme are as follows: a) What is the impact of the self- 
energy feedback and/or the frequency-dependence on the 
type and energy scale of the leading instabilities? Can 
we use the critical scales in this approximation as esti- 
mates for gap sizes in the ordered state at lower scales ? 
b) What can be learned about the flowing self-energy? 
How is its frequency and wavevector-dependence for typ- 
ical situations in the one-band Hubbard model? Is there 
any non-Fermi liquid physics that can be deduced from 
these flows? 



II. MODEL AND METHOD 

Throughout this paper we study the two-dimensional 
Hubbard model on a square lattice with nearest- and 
next-nearest-neighbor hoppings t and t' , respectively, and 
a simple onsite interaction U . The Hamiltonian is given 
by 



= £>(k)4(k)c ff (k) 

k,<7 



E 



(1) 



Here cj.(k) and c a (k) are creation and annihilation op- 
erators with momentum k and spin a € {t, -1} and 
ni,cr = c\ a Ci^ is the density operator at site i. The dis- 
persion relation reads as 

e(k) = — 2i(cos k x + cos k y ) + At' cos k x cos k y — /i , (2) 

where fj, is the chemical potential and the lattice constant 
is set to unity. 

We employ the fRG flow equations for the one-particle- 
irreducible (1PI) vertex function, which are well ex- 
plained in detail in the literature (see, e.g., Refs. fTl f22lj23h . 
The aim is to obtain an effective theory for the low en- 
ergy degrees of freedom of the system at hand. To do 
this, one introduces a regulator or cutoff function 9 A (k) 
in the quadratic part of the fermionic action which de- 
pends on the RG scale A so that the free Green's function 
becomes 



n 



6> A (k) 



(3) 



with the Matsubara frequency u. For the cutoff function 
# A (k) we choose here a momentum-shell cutoff, 



l (k) = 9(| e (k)|-A) 



(4) 



where is the step function. In the numerical imple- 
mentation, the step function is slightly softened for bet- 
ter feasibility. The cutoff function suppresses all modes 



with |e(k)| below the scale A. At the beginning of the 
fRG flow at A = A > max(|e(k)|), i.e. greater than the 
bandwidth, the vertex function of the theory take their 
bare values, as all perturbative corrections suppressed by 
the cutoff function. During the fRG flow the scale A is 
lowered step by step and in principal we recover the com- 
plete action with all vertex corrections at A = 0. 

With the so-introduced dependence of the free Green's 
function on the parameter A one can derive a differen- 
tial equation for the generating functional of the 1PI 
vertices^. After an expansion in the Grassmann fields 
one arrives at an integro-differential equation for all 1PI 
vertex functions. In practice we truncate the infinite hi- 
erarchy of flow equations by setting the three-particle 
(six-point) vertex and all higher vertices to zero. In the 
bare action at A = the three-particle vertex vanishes 
exactly. In principle, during the flow it attains a finite 
value, which is however neglected in our approximation. 
Various fRG studies on the Hubbard-model have shown 
that one can still expect reasonable results as long as one 
restricts the study to weak or moderate coupling. If the 
temperature is low enough, a certain part coupling func- 
tion will typically develop very large values compared to 
the bandwidth before all modes are integrated out, i.e. 
will undergo a flow to strong coupling at a nonzero criti- 
cal scale A c . Here we have to stop the flow as neglecting 
the higher-order vertices is no longer justified. The crit- 
ical scale A c can be used as estimate for ordering tem- 
peratures and concomitant spectral modifications, e.g. 
through gap openings. 

With this approximation and in presence of spin rota- 
tional symmetry the fRG differential equation amounts 
to the following equation for the self- energy £ A and two- 
particle vertex V A , 



— £ A (fc) = jdk'S A {k') V A {k,k',k')~2V A (k,k',k) 



dA 



V A {k ll k 2 , k 3 ) = Tpp + t a H4 + r A HyCr , 



(5) 
(0) 



with the particle-particle channel 

Tp P (k 1 ,k 2 , fc 3 ) = 

- / dkV A {k u k 2 ,k)L A {k,qp P )V A {k,q PP M), (?) 



the direct particle-hole channel 



TpH,d(k 1 ,k 2 ,k 3 ) = - dk 



- 2V A {k ll k, k 3 )L A {k, qPH,d)V A {q P H,d, k 2 , k) 

+ V A (k, k u k 3 )L A (k, qPH,d)V A {q P HA, k 2 , k) 

+ V A (k 1 ,k,k 3 )L A {k 1 qp H , d )V A (k 2 ,q PH ^k)\ , (8) 
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and the crossed particle-hole channel 
T A p H M k i,k2,k 3 )= (9) 

J dkV A (k,k 2 ,k 3 )L A {k,qpH,cr)V A (k 1 ,qpH,cr,k) , 

with the combined index k = (k, ui). Here, qpp = (— k + 
ki+k 2 , -W + W1+W2), qph.d = (k+ki-k 3 ; W + W1-W3), 
qph,cr = (k + k 2 — k3,u> + W2 ~ W3) are the quantum 
numbers of the second loop line. We used the shorthand 
notation for the momentum integral and Matsubara sum- 
mation J dk — (27)7 f dkJ2 u > with the temperature T. 
The fourth frequency and momentum of the interaction 
vertex as well as of the second internal line are fixed by 
conservation. The spin convention is such that the spin 
is the same for the first and third leg. The internal loop 
is given by 

L A (fc, k') = S A (k)G A (k') + G A (k)S A (k'), (10) 

with the full propagator G A (k) = [(G^(fc))- 1 - S A (/c)]- 1 
and the single-scale propagator 

S A (k) = -G A (k) (^[Gft*)]- 1 ) G A (k) 

^~ e (k)]d A fl A (k) 

[iw - e(k) - A (k)£ A (fc)] 2 ' 1 ' 

Note that we do not employ the Katanin- modification 24 , 
i.e. the replacement of S A (k) by ■j K G A {k) in Eq. (JTDJ, 
as this would increase the numerical effort very strongly 
in this two-dimensional problem. 

At the beginning of the flow the coupling function as- 
sumes the value of the onsite interaction V A {k\, k 2 , k 3 ) = 
U and the self-energy vanishes E A (fc) = 0. 

Here we want to study the fRG flow with frequency- 
dependent self-energy and two- particle vertex. To do 
this, we divide the Brillouin zone into Nk patches to 
cover the momentum-dependence of the vertices at the 
Fermi surface according to Ref. [3 (see Fig. CD). This 
amounts to a projection of the momentum-dependence 
onto the Fermi surface and the vertex has then three 
momentum indices k^ = 1, ...JVfe. Additionally, we di- 
vide the Matsubara axis in N u patches and calculate the 
Matsubara sum numerically. In this work, we use two 
different patching schemes. In the first scheme (Dl) the 
vertices and self-energy are calculated at the 10 Matsub- 
ara frequencies with the smallest absolute values. The 
Matsubara frequencies of the respective patches are thus 
given by u p = {±1, ±3, ±5, ±7, ±9} x ttT. We keep the 
smallest frequencies as theses are expected to give the 
leading contributions to the flow. Additionally we con- 
sider a patching where the frequencies of the patches are 
at cup = {±1, ±3, ±5, ±11, ±21} x ttT to cover the de- 
pendence at higher frequencies more accurately. This 
discretization is referred to as D2. In both schemes the 
vertices and self-energy at all other frequencies occur- 
ring in the loops are assumed to have the same value 
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FIG. 1: N-patch discretization scheme with — 32. The 
Fermi surface with t' = 0.2t, fi = —At' is plotted in red. The 
patches are divided by black dotted lines. The red dots denote 
the momenta at which the coupling function and self-energy is 
computed. The umklapp-surface is shown as the blue dashed 
line. 



as the closest frequency-patch. This amounts to keeping 
the frequency-dependence in self-energy and interactions 
constant above the largest positive or smallest negative 
discretization frequency. 

In the loop diagrams with self-energy inclusion, the 
Matsubara sums have to be perfumed numerically. To do 
this, we truncate the Matsubara sums at some Woo and 
check whether further increasing this frequency cutoff 
changes our results. In practice we chose = 15007rT, 
as further increasing the maximum frequency did not lead 
to major changes in self-energy or critical scales. 

This way we calculate the vertex and self-energy only 
for a given set of Matsubara frequencies and momentum 
vectors on the Fermi surface. The flowing self-energy 
then appears on the internal lines of the loop diagrams, 
as visible in Eqs. (ITU|) and (fTTj) . Importantly, we neglect 
the real part of the self-energy. This is because a feed- 
back of the real part of the self-energy would lead to a 
Fermi surface shift which is hard to handle within the N- 
patch scheme, as the Fermi surface could be shifted into 
a momentum region which has already been integrated 
out in previous fRG steps, causing severe singularities in 
the flow. 



III. RESULTS 

A. Case without instability 

First, as a sanity check in order to get familiar with the 
data we obtain from the fRG, let us look at a case where 
the flow indicates normal metallic behavior. For that we 
choose t' = 0.25t and /j, = —2t such that the Fermi sur- 
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face is away from the van Hove situation and from the 
antiferromagnetic (AF) spin-density wave (SDW) insta- 
bility. Then, if we stay away from lowest temperatures, 
the fRG flow does not lead to strong coupling, i.e. we can 
integrate out all modes without encountering a strong 
growth of any components of the coupling function. In 
Fig. [2] we show the imaginary parts of the self-energy as 
function of the Matsubara index for different tempera- 
tures. Besides a decrease of the self-energy toward lower 
T, we observe a monotonic frequency-dependence that 
features some kind of downward step around zero fre- 
quency The slope of the overall curve can easily be in- 
terpreted as quasiparticle weight or Z-factor, while the 
step can be understood as a finite lifetime r. 

Note that for higher frequencies of the scale of the 
bandwidth, the imaginary part of the self-energy should 
approach zero again, but our frequency window is not 
large enough to capture this behavior. In our treatment, 
the self-energy at frequencies outside the discretization 
window (i.e. for Matsubara indices higher than ±10) is 
approximated by the value at the closest discretization 
frequency. This may underestimate the true self-energy 
in the frequency range where the linear slope around 
zero frequency still continues, but may overestimate it 
at higher frequencies where the imaginary part of the 
self-energy is a decreasing function of the absolute value 
of the frequency. To some extent these two effects will 
cancel each other, and the large-frequency self energy 
is anyway less relevant for the Fermi surface instabili- 
ties considered here. Therefore we think that our simple 
scheme of treating the self-energy behavior at higher fre- 
quencies is qualitatively adequate. 

We can analyze the behavior of the quasiparticle life- 
time obtained from the step around zero frequency a bit 
more closely. In a Fermi liquid, the life-time is related 
to the temperature as r _1 ~ T 2 , with a logarithmic cor- 
rection in two dimensions that we might not be able to 
resolve. Let us test if the law 



ocT" 



(12) 



is verified by our data, where we should find a = 2 for 
a Fermi liquid with round Fermi surface. At van Hove 
filling due to the divergence of the density of states the 
exponent is expected to be altered to a = 1—. Note 
however that all these literature values are obtained by 
finite-order perturbation theory, while the fRG sums up 
infinite orders in the bare coupling. The fRG estimate 
for the life-time is obtained by linearly extrapolating the 
self-energy for the two lowest frequencies to uj = 



= AE A = ^Lm[E A (k F ,3-KT)} - ^Im[S A i 



The expression 

AClogCr- 1 )) 



k F ,wT)) 
(13) 



log^ 1 ) 



A(log(T)) log(T s 



i+ll 



log(li) 



(14) 



gives an estimate for the exponent, extracted from data 
at two temperatures and Tj+i. The right plot of Fig. 
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FIG. 2: Left: Imaginary part of the self-energy in units of 
t for different T > T c . t' = 0.25t, p = -It. Shown are 
the values at momentum k = 1 (see Fig. [T]). The largest 
couplings at A = are of the order of 4t. Right: Estimation 
for the exponent in Eq. (|12p for t' = 0.25 at van Hove filling 
(red triangles) and away from van Hove filling (blue squares, 
green diamonds), obtained via Eqs. (| 13|) and (|14|) for k = 1 
(see Fig. [1} . All data for U = 3t with self-energy feedback in 
the discretization Dl. 



[5] shows the temperature-dependence of a. Although we 
do not observe a constant a presumably due to the error 
of the extrapolation, at least a for the rounded Fermi 
surface is near the expected value of two for low enough 
temperatures, and similarly for van Hove filling a is close 
to one for low temperatures. 



B. Critical scales and leading instabilities 

Now let us come to cases with flows to strong cou- 
pling. The first question we want to ask is whether the 
inclusion of frequency-dependence and feedback from the 
imaginary part of the self-energy leads to a different lead- 
ing instability or a change in the critical scale. Here we 
study this question along the parameter line of van Hove 
fillings where the Fermi surface contains the points (ir, 0) 
and (0, 7r). These are saddle points of the dispersion and 
cause a divergence of the density of states. From previ- 
ous works it is known that this van Hove situation leads 
to competing ordering tendencies. Hence this situation 
is a good arena to study the impact of the inclusion of 
frequency-dependence and self-energy feedback. 

First let us look at the van Hove situation at half filling, 
where in addition the Fermi surface is perfectly nested. 
Previous fRG works have clearly established an AFM in- 
stability, signaling an AF SDW ground state for this case. 
This in in accordance with the expectations from simple 
RPA arguments, but also with quantum Monte Carlo^. 

For t' = 0, /J, = 0, the character of the instability in 
the fRG does not change when we include a frequency- 
dependence of the effective interactions, or when we in- 
clude the self-energy feedback. This means that the 
couplings that drive the antiferromagnetic susceptibility 
with wave vector transfer (tt, tt) diverge most strongly, as 
described in earlier worksi. In Fig. [3]we show the critical 
scales obtained in the fRG in the various approximations 
(without frequency-dependence of the effective interac- 
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FIG. 3: Critical scale, at which the couplings exceed 20i, vs. 
U for the fully nested case, with T = O.Oli, t' = 0, and [i — 
for the flows with frequency independent vertices without self- 
energy feedback (blue squares), with frequency-dependent 
vertices and without self-energy feedback (red triangles), with 
frequency-dependent vertices and with self-energy feedback 
(green diamonds). 



tions and without self-energy feedback, with frequency- 
dependence of the effective interactions but without self- 
energy feedback, with frequency-dependence of the effec- 
tive interactions and with self-energy feedback) as func- 
tion of the interaction strength U. These scales do not 
differ much quantitatively. One way to understand the 
agreement is that the instability in this case is very clue 
to a standard RPA instability, where only the most sin- 
gular frequency transfer matters, and different frequency 
transfers do not couple, hence the frequency resolution 
does not come into play. Furthermore, the flow goes off to 
strong coupling at scales before the self-energy becomes 
noticeable. Neverthelss, with respect to plain RPA, the 
critical scales are reduced by the channel coupling in the 
fRG. Hence, in this case, the simplest fRG approximation 
without frequency-dependence and without self-energy 
feedback is already quite good and cannot be improved 
much within the range of numerical possibilities. The 
remaining approximations are besides the discretizations 
of wave vector and frequency-dependences the neglect of 
the real part of the self-energy and the truncation of the 
flow hierarchy after the four-point function. We expect 
that the any remaining differences are due to these two 
approximations. In particular for larger U ~ 4t, there 
should be precursors of the spectral weight transfer that 
ultimately leads to the opening of a Mott gap at U ~ 8t. 
As is believed that the truncated fRG does not allow the 
description on the Mott transition, it is plausible that the 
spectral weight transfer is still not described correctly in 
our improved approximation. Furthermore, the physics 
of collective fluctuations might not be captured to some 
extent due to the truncation after the four-point vertex. 
Hence, if we compare roughly to gap scales found in non- 
perturbative cluster calculations for large clusters^, the 
A c -curve seems to rise too steeply as function of U. How- 




FIG. 4: Critical scale A c , at which the couplings exceed 20t, 
vs. next-nearest-neighbor hopping i'/i at van Hove filling for 
U — 3i with N u = 1 and without self- energy feedback (blue 
squares), with frequency-dependent vertices and without self- 
energy feedback (red), with frequency-dependent vertices and 
with self-energy feedback (green). Two different frequency 
discretization Dl (triangles, solid lines) and D2 (diamonds, 
dashed line) are used. Left: T = O.Oli. Right: T = A°' 01 is 
the the critical scale for T = O.Oli without self-energy feed- 
back with discretization Dl (red solid line in the left plot). 



ever, regarding the critical scale as estimate of the gap 
scale, up to a factor of order one, the fRG is in the right 
range. 

Let us now look at the situation where the particle- 
hole nesting is destroyed by larger t'. In FigHJ we plot 
the critical scale vs. the next-nearest-neighbor hopping t' 
at van Hove filling (i.e. chemical potential fi = —At'), at 
which the van Hove energy with a logarithmic diverging 
density of states is at the Fermi level. The fRG flows 
in Fig [3] a) have been computed at T = O.Oli. Lowering 
the temperature further increases the numerical effort, as 
the maximum cutoff frequency uj x then has to be taken 
larger too. 

As can be seen in Fig. 2] a) for T = 0.014 the inclusion 
of frequency-dependence and the different discretizations 
affects the critical scale only slightly. The self-energy 
feedback however reduces the critical scale significantly 
especially for larger t'. We also compared how differ- 
ent patchings on the Matsubara axis affects the critical 
scales. In Fig. 2]we show data for two different discretiza- 
tions of the Matsubara axis. At low temperatures, the 
two different discretization lead to very similar results, 
while at higher temperatures the two discretizations show 
some quantitative differences, depending on t' . 

One might suspect that the critical scale for higher 
t' is more strongly suppressed because the critical scale 
comes closer to the chosen temperature, so that the dif- 
ference in the critical scale is merely a finite temperature 
effect. To check this we plot the critical scale for differ- 
ent t' with T = A° 01 , where A° 01 is the critical scale at 
T = O.Oli without self-energy feedback. Now the ratio 
of temperature and critical scale is comparable for all t'. 
We see that inclusion of self-energy-feedback suppresses 
the critical scale now already for small t' , but definitely 
more strongly at larger values of t'. 

Interestingly we find that increasing the temperature 
reduces the critical scale more strongly, when frequency- 
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dependence is neglected at small t', so that we discover a 
parameter space, in which now the critical scale increases 
when one includes frequency-dependence in contrast to 
Fig. 2]a). At larger temperatures the difference between 
the two discretization schemes also becomes more pro- 
nounced. The flows with the discretization D2 now ex- 
hibit a larger critical scale than the flows with Dl. Hence 
the precision of the description is better at low T. 

We want to have a closer look at the origin of the 
behavior observed in Fif. 2J A natural question to 
ask is whether the AFM SDW and the superconduct- 
ing (SC) tendencies respond differently to the inclusion 
of frequency-dependence and self-energy, and whether 
the tentative phase diagrams drawn by determining the 
leading sintability are changed by these improvements of 
the scheme. We therefore consider the effective coupling 
strengths of the AFM and SC channel, which we define 
as 



Xafm = — yA KP' (k = 1, wo), (k = 24, uj Q )) , 
k p 

(15) 

*SC - ^2 E ^ « k ' W °)> (~ k > " W o), (P> /( k )/(P) 
k,p 

(16) 

Here, k = 1 and k = 24 denote the patches with momen- 
tum vectors closest to the van Hove singularities roughly 
connected by momentum transfer (jr, 7r) (see Fig. [TJ and 
ujq = ttT. /(k) and /(p) are <i-wave form factors and 
are given by /(k) = cos(k x ) — cos(k y ). We consider the 
vertices of the smallest Matsubara frequencies, as these 
grow most strongly during the fRG flow. If these aver- 
ages diverge, the corresponding susceptibilities diverge as 
well, so that these quantities can be used as a measure of 
the coupling strengths of the respective channel. Typical 
flows of these quantities are shown in the left plot of Fig. 
[5j As a criterium for the leading instability, we choose 
the derivatives of the coupling strength of the effective 
channels with respect to A at A c . A c is again defined as 
the scale at which the couplings exceed a value of 20t. 
The results are shown in the right plot of Fig. [5j The 
intersection of the respective curves for the AFM and SC 
coupling strength is taken as an estimate for the bound- 
ary between the regimes. Of course the precise values 
from this procedure depend somewhat on the chosen def- 
inition of A c . On the other hand, this procedure suffices 
for the more qualitative discussion, which we are inter- 
ested in, namely the impact of frequency-dependence and 
self-energy effects on the competition between the differ- 
ent orders. The comparison for the different cases, i.e. 
with and without frequency-dependence and self-energy 
feedback, is summarized in table HI 

At T — O.Olt and without self-energy the bound- 
ary between the regimes lies at t' ss 0.1 9t and without 
frequency-dependent vertices. If we include the feedback 




FIG. 5: Left: Flow of the effective strength of the antifer- 
romagnetic (solid lines) and superconducting (dashed lines) 
channel given in Eq. (|15l I16|) in units of t for t' = 0.2t 
(blue) and t' = 0.26t (green) in the discretization scheme Dl 
and with self-energy feedback. Right: Slope of the effective 
strength of the antiferromagnetic (solid lines) and supercon- 
ducting (dashed lines) channel at the scale A c , at which the 
couplings exceed 20t. Shown are the data with frequency in- 
dependent vertices and without self- energy feedback (blue 
squares), without self-energy feedback (red triangles), and 
with self-energy feedback (green triangles), both in the dis- 
cretization scheme Dl. The intersection point of the two 
curves is used as an estimate for the boundary between the 
AFM and SC dominated regime. All data with U = 3t, at 
van Hove filling and T = O.Olf. 
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TABLE I: Approximate values for t'/t, for the transition from 
the AFM to the SC regime at van Hove filling and U = 3i. 
The values are obtained by the comparison of the slope of the 
effective strength of the respective channels at the scale A c , 
at which the couplings exceed 20t (see also Fig. [5}. 

of self-energy this boundary is shifted towards higher t' . 
We conclude that life-time effects harm the supercon- 
ducting instability more than the AFM instability. 

For the second run at higher temperatures T = A° 01 
the phase boundary is generally shifted towards higher 
t', too. This basically reflects the fact that the relevant 
energy scale of the superconducting instability is lower 
than the one of the AFM instability. Hence the pair- 
ing channel is more strongly affected by the finite tem- 
perature. The difference between the different levels of 
approximations is much larger than at T = O.Olf in cor- 
respondence with the larger discrepancy in the critical 
scales in Fig. Furthermore, we note that the differ- 
ent discretization schemes yield qualitatively similar re- 
sults, which confirms that the picture drawn here for the 
leading instabilities is fortunately rather independent on 
details of the numerical implementation. 

The overall conclusion we can draw from this analysis 
is that the frequency-dependence of the couplings and 
the self-energy effects are certainly important for quanti- 
tative questions. However, they in general do not qualita- 
tively change the physical picture obtained from the sim- 
pler flow schemes with frequency-independent couplings 
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FIG. 6: Imaginary part of the self-energy in units of 4 for T > 
T c at A = for three parameter choices at van Hove filling, 
U = 3t. T — OM, t' = (blue); T = OM, t' = 0.254 (red); 
T — O.lt, t' — 0.254 (green). Shown are different patches, 
i.e. momenta at the Fermi surface k = 1 (solid line), k = 2 
(dashed line), k = 3 (dashed dotted line), k = 4 (dotted line), 
all nearly on top of each other: only weak angular dependence. 
The left and right plot shows the same data with different 
scaling of the frequency axis. All data with discretization 
scheme Dl and feedback of self-energy. The largest couplings 
at A = are of the order of 104 for T = OM, t' = 0.254 (red). 
The parameter choices shown in blue and green are close to 
the critical temperature, here the largest couplings are close 
to 204. 

and with neglected self-energy. 



C. Flow of the self-energy 

Next we consider the flow of the self-energy for finite 
temperatures above T c . The main goal is now to analyze 
the frequency- and wavevector-dependence. 

We run the fRG flow for different T > T c such that 
all modes can be integrated out. Practically, the flow is 
stopped at a scale A = 2 x 10~ 4 i, which is much lower 
than the temperature. Below these scales, the vertices 
and self-energies do not get renormalized substantially 
by lowering A even further. We stay in the weakly to 
moderately coupled regime at all scales. 

In Fig. [5] the imaginary part of self-energy is plotted 
over the Matsubara frequency for different momentum 
patches, for the fully nested case with t' = at half fill- 
ing and for van Hove filling with t' = 0.254. There is 
only a weak angular dependence at these elevated tem- 
peratures above the instability, independent of the Fermi 
surface shape. We checked that the weak angular depen- 
dence is reproduced in the other frequency discretization 
schemes as well. On the left hand side of the plot, we 
show the data as function of the Matsubara frequency, 
i.e. with different frequency window oc T for the two 
temperatures. On the right hand side we show the same 
data vs. iui/(nT) as function of the Matsubara frequency 
index. 

At higher temperatures T — 0.54, we observe again 
a step-like discontinuity in the self-energy at zero fre- 
quency on the Matsubara axis. This can be interpreted 
as an inverse lifetime of the quasiparticle peak. If we go 
to lower temperatures but remain above the instability 
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FIG. 7: Imaginary part of the self-energy in units of t for 
different discretization schemes Dl, D2 at the scale Aiot at 
which the couplings exceed a value of lOt. Aiot — 0.0844 
(Aiot = 0.0934) for Dl (D2). Shown are different momenta 
at the Fermi surface k = 1 and k = 4 (see Fig. [TJ. All data 
for T = 0.014, t' = 0.24, U = 3, at van Hove filling and with 
feedback of self-energy. 



(now only for the curved Fermi surface, as otherwise the 
critical scale is too high), this step gets much smaller 
(also because the Matsubara frequencies move close to- 
gether) and we are left with a rather linear frequency- 
dependence, which can again be captured by a quasipar- 
ticle weight Z. This again turns out to be rather inde- 
pendent of the location on the Fermi surface. So all this 
data at fixed temperature looks quite consistent with a 
rather normal and isotropic metallic system, at least at 
these temperatures above the instability. 

There have been a number of more detailed studies 
(e.g. Refs. l28U29f ) of the self-energy in the case of van 
Hove filling, pointing out special, non-Fermi liquid-like 
properties due to the diverging density of states. At a 
given T, it is difficult to read out any particular self- 
energy property from our data at t' = 0.254. However, 
if we measure the temperature-dependence of the step- 
like feature that we interpreted as inverse lifetime and 
fit it to the law T a , we obtain a different exponent with 
a approaching 1 toward lower T. This marginal Fermi- 
liquid behavior-2S is shown in Fig. [5] and consistent with 
second-order predictions for the van-Hove situatio n 25 i 31 . 

In Fig. [7] we show data lower temperatures for the 
frequency-dependence of the self-energy, now comparing 
different frequency discretizations. We compare data at 
a scale where the couplings exceed a value 10i. Note that 
this is at slightly different scales in both cases. We see 
that for low frequencies the discretizations compare quite 
well. At larger frequencies the quantitative agreement 
becomes worse. However, the qualitative features, as for 
instance a larger slope of the self-energy near the saddle 
points compared to the Brillouin zone diagonal, is found 
consistently in both schemes. Quite generally, our data 
shows that the self-energies in the two schemes are more 
consistent at lower temperatures. 
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For T <T C the flows go to strong coupling. If we stop 
the flow when the vertices exceed a value of the band- 
width, the momentum-dependence of the self-energy are 
rather weak. However, if we continue to flow towards 
the instability, stronger anisotropics are found. Strictly 
speaking, for these large values of the vertices the trunca- 
tion of the fRG flow equation with neglect of the higher- 
order vertices is no longer justified. We expect however 
that we can still get some qualitative insight in the break- 
down of the Fermi liquid behavior for low temperatures, 
similar to the argumentation in previous work s 12 ' 13 . The 
corresponding data is shown in Fig. [SJ both with and 
without self-energy feedback. 

In general, the frequency-dependence of the imaginary 
part of the self-energy at these low temperatures near the 
instability is dominated by a linear decrease that gets 
steeper the closer we get to the critical scale. We can 
measure this angle-dependent slope £' A on the frequency 
window and encode it in a quasiparticle weight 



Z A (k F ) 



[1 - S' A (k F ) 



where we define E' A as 

E A (k F ,cj m ) - S A (k_F, -uj r , 



S' A (k F ) 



2uj, 



(17) 



(18) 



Here E A (k F ,±u; m ) are the self-energies at the Matsub- 
ara frequencies with largest absolute value uj rn — 9ttT 
available in discretization Dl. The flow of these weights 
is shown in Fig. [9l We see that the if-factor diminishes 
toward the instability, but in a rather weak, almost loga- 
rithmic way, such that this is possibly not the most severe 
effect in the self-energy. Besides this slope, we also notice 
again a step AE A around zero frequency (Eq. [13]) , that 
grows rapidly toward the instability, as shown in the right 
plot of Fig. [SJ We also show curves for different locations 
on the Fermi surface. These different lines indicate the 
growth of the anisotropy toward the instability for the 
van Hove situation away from half filling, with stronger 
self-energy effects for Fermi surface points near the van 
Hove points. In the perfectly nested case, the anisotropy 
is less pronounced. 

Previous works on the flow of the self-energy using dif- 
ferent approximations have arrived at some similar but 
also some different conclusions regarding the self-energy 
flow. In Refs. UH ancU 3 - the anisotropy appears to be 
somewhat more pronounced in the sense that the Fermi 
surface point closest to the van Hove points showed in- 
dications for a break-up of the quasiparticle peak into 
two peaks, with a spectral gap opening between them. 
For this to occur, one needs an additional low-energy 
feature in the self-energy. On the imaginary axis, this 
shows up as an increase of the absolute value of the 
self-energy imaginary part towards low frequencies, that 
turns around the linear slop and eventually leads to a 
l/(«w)-pole. In Ref. [lathis could only be seen clearly by 
interpolating between the lowest Matsubara frequencies, 
an option we cannot choose in our flow with frequency- 
dependent coupling. However, the growing step AE A can 
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FIG. 8: Imaginary part of the self-energy in units of t for 
t' = (blue) and t' = 0.2 (red) at van Hove filling for U = 3t 
at the scale A&ot where the couplings exceed a value of 80t. 
A 80 t = 0.26t (A 80t = 0.019t) for t' = (t' = 0.2t). Shown 
are different momenta at the Fermi surface k = 1 (solid line), 
k = 2 (dashed line), k = 3 (dashed dotted line), k = 4 (dotted 
line) (see Fig. [TJ. All data with discretization scheme Dl and 
with feedback of self-energy. 



be interpreted as a precursor of this effect, as it leads to a 
breaking up of the simple linear decrease. It is well pos- 
sible that these subtle effects very close to the instability 
come out in different extent in the different approxima- 
tions. Also in Ref. [H, where the self-energy is evaluated 
directly on the real frequency axis, a new low-energy fea- 
ture in the self-energy emerges very close to the insta- 
bility and eventual splits the quasiparticle peak in two. 
Another aspect of this feature is a rise in the scattering 
rate, which is consistent with our increasing step AE A . 
On the real axis, the linear slope of of the self-energy 
real part goes into the Z-factor. In Ref. [l3|, this slope 
after subtracting the additional low-energy feature is not 
strongly scale-dependent, very much like our data, where 
the Z-factor do not really dive to zero toward the insta- 
bility. Since the precision of the truncated flows is not 
very high in this regime, we refrain from a more elab- 
orate discussion of these differences. One message one 
may learn from these studies is however that the self- 
energy in this regime should not only be parameterized 
with a Z-factor. At least close to the instability, it devel- 
ops additional low-energy structures besides this linear 
envelope captured by the Z-factor. 



IV. CONCLUSIONS 

We have analyzed the low-energy properties of effec- 
tive interactions and self-energies in the two-dimensional 
Hubbard model using a fRG approach. In contrast with 
previous approaches to the same model ( e.g. Refs. 
and many more works cited in Ref. [l]), we take into ac- 
count the frequency-dependence of the interactions and 
the feedback of important parts of the self-energy on the 
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A/t 



FIG. 9: Flow of the Z-factor as given by Eqs. JT71 [18]) and 
the step AE A of the self-energy in units of t at oj — defined 
as in Eq. (|13p for different momenta at the Fermi surface (see 
Fig. [[]). U = 3t, t' = 0.2t, van Hove filling, T = O.Olt, with 
discretization scheme Dl and with feedback of self-energy. 



flow of the latter. For this, we discretized the Matsubara 
frequency axis in different schemes with N u = 10 patches. 
For the self-energy we took into account the imaginary 
part (on the Matsubara axis) only, while we ignored the 
real part. The argument for ignoring the latter is mainly 
feasibility, a flowing dispersion is technically difficult to 
deal with, and would require additional devices to keep 
the particle number fixed. Our experience from previous 
studies^ is that the the impact of the dispersion renor- 
malization on the instability is mainly quantitative and 
does not alter the regime found from the flows to strong 
coupling. 

The main goal of this study was to see how the fRG 
flows of the simpler studies are altered by the straightfor- 
ward inclusion of frequency-dependence and (imaginary 
part of the) self-energy feedback. Here, after a first san- 
ity check on a case without flow to strong coupling, we 
analyzed how the critical scales for the flow to strong 
coupling change, and how the character of the instability 
in terms of the leading ordering tendency changes de- 
pending on the different approximations. The overall re- 
sult is that the frequency-dependence and the imaginary 
part of self-energy usually decrease the critical scales by 
a factor at most of order one. Only near quantum crit- 
ical points, where the critical scales vary strongly and 
become small, the change can be more drastic. Further- 
more, the competition for the leading instability is also 
affected only quantitatively. We find the same structure 
of the phase diagrams and have not seen that any of the 
regimes found previously are wiped out by these alter- 
ations of the scheme. In particular there is a sizable 
regime with dominating d-wave pairing in the ground 
state. 

Above the instability temperatures, the self-energy 
did not show strong signatures in its frequency- or 
wavevector-dependence. The basic structure could be 
captured by an inverse lifetime or scattering rate and a 
Z-factor for the quasiparticle weight. These parameters 
turned out to be weakly wavevector-dependent around 
the Fermi surface unless one gets close to the instabil- 
ity at lower temperatures. At van Hove filling and away 
from perfect nesting, the scattering rate above the insta- 



bility comes out consistent with marginal Fermi liquid 
behavior, as found in previous studies^!. This how- 
ever does not prevent the low-T instabilities from occur- 
ring, i.e. even if the self-energy behaves different from a 
Fermi liquid, the instabilities found without inclusion of 
the self-energy are still present. 

These findings on the behavior of the self-energy above 
the instability make a lot of sense, but do not reveal any 
unexpected physics. Certainly, the self-energy flows cor- 
roborate second-order -perturbation theory results on the 
non- Fermi-liquid behavior at the van Hove filling, and the 
low-T instabilities are not changed qualitatively. Hence 
the present work represents an important check for the 
body of knowledge has been obtained using simple fRG 
flows without self-energy feedback (for a review, see Ref. 
1). Basically, the message we infer from the data is that 
frequency-dependence and self-energy effects are neces- 
sary to deal with if one is trying to obtain truly quanti- 
tative results, but on the qualitative level, regarding the 
leading instabilities and orders of magnitude, these im- 
provements do not lead to significant changes. So the 
simple flows without self-energy are expected to be good 
guides through the basic phase diagrams. Similar state- 
ments also hold for parameter trends, e.g. as function of 
the chemical potential of the band structure parameters. 
Here, we have found that these trends occur irrespective 
of the approximation level. 

In this paper, we do not discuss explicitly the fre- 
quency structure of the effective interactions. This is 
quite rich and contains some interesting aspects, some 
of which has been discussed without self-energy feedback 
in Ref. [l9|. Here, it is however difficult to describe the 
data with a few simple parameters, and less literature 
is available for comparison. Hence we did not elaborate 
on this here. It would however be interesting to com- 
pare the frequency-resolved vertex functions from this 
weak-coupling approach with data obtained with QMC 
or dynamical mean-field-like cluster methods, in order 
to assess the effects of strong correlations (possibly diffi- 
cult for the fRG) and larger distances (possibly less well 
captured by the strong-coupling approaches) on these ob- 
jects. 

Regarding the theoretical approach chosen here, we 
were able to show that it leads to useful results. On the 
other hand, in order to look for subtle effects near the in- 
stability such as pseudogaps or more detailed non-Fermi 
liquid behavior, the chosen approach appears somewhat 
to heavy and inflexible. The numerical effort due to a 
vertex that depends on three wavevectors and frequen- 
cies makes it rather cumbersome to compare different 
discretizations or implementations. Here, some clever re- 
duction of the information carried along my serve useful. 
Our current code my serve as a guide how good these 
alternatives are. 

At the time of writing of this paper, another work 
treating the self-energy feedback in fRG flows in the two- 
dimensional Hubbard model has appeared as a preprint 33 
which may point in the right direction. This study em- 
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ploys the channel decomposition of the interaction vertex 
proposed by Husemann and Salmhofer™. The restriction 
to a few form factors of the interaction channels then 
allows the authors to reach a higher momentum space 
and frequency precision as in our straightforward patch- 
ing scheme, and to safely go to lowest temperatures and 
frequencies. Furthermore they also consider the flow of 
the real part of the self-energy which was ignored in our 
scheme for simplicity. The choice of a soft frequency cut- 
off instead of a momentum-shell cutoff picks up ferromag- 
netic tendencies, which are ignored in our approach and 
which would be important for t' > 0.3i near van Hove 
filling. So the results of both schemes should only be 
compared at small t' . There seems to be good agreement 
at least in two aspects, which are the suppression of the 
critical scales due to self-energy effects, and the stability 
of the d-wave pairing regime with respect to inclusion of 
these aspects. Ref. [2l| also describes a non-Fermi liquid 
behavior of the self-energy at the van Hove-filling with 
exponent 3/4 at t' = 0.355i, which is slightly less than 
our second-order-like exponent near one determined for 
t' = 0.25f. In comparing the pros and cons of these two 
related schemes, it is important to notice that our direct 
patching implementation without channel decomposition 



does allow the vertex to depend on three different wave 
vectors and Matsubara frequencies, while the approach 
in Ref. [33| only allows for certain form factors in momen- 
tum space and treats the frequency-dependence as the 
sum of three functions each depending on one particular 
frequency combinations. However, it is hoped that these 
constraints do not influence the results too strongly. In- 
deed our study shows that taking along more frequency- 
and momentum-structure does not change the main re- 
sults and does not reveal any additional physics. Hence, 
it is quite likely that the channel-resolved formalism of 
Ref. [33| turns out to be a better compromise as it re- 
duces the amount of information carried along in the flow 
in a physical meaningful way. This way one might ob- 
tain more precise data and one might be able to go on 
to more invalid systems without reducing the precision 
again. In a way, the approach chosen in this paper can 
be viewed as complementary, as it is less versatile but 
on the other hand puts less bias in the momentum and 
frequency-dependences of the vertex functions. 

We acknowledge discussions with M. Kinza, A. 
Katanin, K.U. Giering, C. Husemann, M. Salmhofer, G. 
Klingschat and A. Eberlein. This work was supported by 
the DFG research units FOR 723 and 1162. 
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